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We introduce the idea of using adiabatic rotation to generate superpositions of a large class 
of quantum states. For quantum computing this is an interesting alternative to the well-studied 
"straight line" adiabatic evolution. In ways that complement recent results, we show how to effi- 
ciently prepare three types of states: Kitaev's toric code state, the cluster state of the measurement- 
based computation model and the history state used in the adiabatic simulation of quantum circuit. 
We also show that the method, when adapted for quantum search, provides quadratic speedup as 
other optimal methods do with the advantages that the problem Hamiltonian is time-independent 
and that the energy gap above the ground state is strictly nondecreasing with time. Likewise the 
method can be used for optimization as an alternative to the algorithm of Farhi et al pj|. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

The presence of superpositions of classical states lies at 
the heart of many non-classical behaviors of quantum sys- 
tems. It is therefore not surprising that uniform super- 
position states serve as powerful resources for quantum 
information processing. The main result in this article 
is that there is a conceptually simple way of preparing 
a class of linear superpositions by adiabatically rotating 
a driving Hamiltonian. Seen another way, it is also a 
method for finding the ground states of certain Hamilto- 
nians. While the applications we discuss are all related to 
quantum computing, the idea is relevant to preparation 
of quantum states in general. 

Adiabatic evolution studied in quantum computing 
often takes the form of a linear interpolation H(s) = 
(1 -s)H initia i + sHfi na i whether it is for problem solv- 
ing [2 Q or state preparation 0, 0, [Bj]. In this arti- 
cle we consider a different paradigm that is inspired by 
observations in Q and follows a "rotation" instead. A 
time-dependent similarity transform on the entire Hamil- 
tonian, as pointed out in often requires highly non- 
local interactions and thus seems to have few applica- 
tions. One may however adopt a different perspective 
and consider this time-dependent transformation as a 
driving term in the presence of a well-gapped and time- 
independent Hamiltonian. What we will show in Section 
II is that some of the desirable properties associated with 
the similarity transform remains in this more general set- 
ting. 

In Section III, we outline how to use this idea to pre- 
pare states with interesting physical properties. The first 
example is the ground states of Kitaev's toric code 0], 
which have been of much interest because of their nat- 
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ural fault tolerance as memory and exhibition of topo- 
logical order Q. In 31 it is shown that a ground state of 
the toric code Hamiltonian can be prepared adiabatically 
through an linear interpolation in time 0(^/n) where n 
is the number of sites. Both this method and the origi- 
nal preparation-by-measurement method Q are optimal 
since they saturate the Lieb- Robinson boundfiol]. which 
places a theoretical limit on the efficiency of any prepara- 
tion method for topologically ordered states. However it 
requires duality mapping to an Ising model, which means 
the same result cannot be easily reproduced for a differ- 
ent system with topological order. We will show, without 
requiring extra types of interactions, that there is an el- 
ementary, circuit-like adiabatic evolution path that also 
prepares the toric code state optimally. The same idea 
would not only be applicable to other topologically or- 
dered states, but also, for instance the cluster state [l3| 
used in measurement-based computing. Recently it was 
shown in [3| that the a state similar to the cluster state 
in computation power can be prepared efficiently via adi- 
abatic linear interpolation. Our method shows how the 
cluster state itself can be prepared adiabatically and ef- 
ficiently. 

In Section IV, we move on to a more complicated ex- 
ample. In 6] it is noted that the adiabatic linear inter- 
polation corresponding to a quantum circuit used by [B[ 
can be thought of as a special case in a larger family of 
paths. It may still appear, however, that fairly elaborate 
gap analysis [!, |Tl[ would be required to check the effi- 
ciency of any path that generates the history state. This 
turns out to be not the case. We will describe two paths 
that generate the history state efficiently, and apart from 
being amenable to gap analysis, one only requires time- 
dependent interactions on a few qubits. 

In Section V, we turn our attention to problem solv- 
ing. Adiabatic quantum computing was originally pro- 
posed not for state preparation but to solve optimization 
(especially NP-complete) problems. While there is still 
controversy on the value of this approach, it is well-known 
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that the the closely related adiabatic quantum search 
method [2| yields an optimal result, in the sense that 
it is as efficient as Grover's algorithm. Here we present 
an adiabatic rotation version of quantum search. Like 
the linear interpolation version, it is also optimal. The 
difference is that in our case, the spectral gap between 
the ground state and the first excited state is strictly 
nondecreasing, and this feature can be extended to other 
optimization problems. Clearly this would make it eas- 
ier for us to determine whether a given problem can be 
solved efficiently. Our approach, however, has an addi- 
tional difficulty compared to linear interpolation. While 
exact examples of linear interpolation algorithms often 
assume nonlocal interactions (as [2] does), one can al- 
most always formulate a linear interpolation algorithm 
with local interactions. The situation is different with 
adiabatic rotation, because it requires separation of the 
initial state from the computational states. We will con- 
clude with a brief discussion of what this means and some 
suggestions for future direction. 



II. THE BASIC FORMULATION 

In the simplest case, our Hamiltonian consists of a 
time-independent piece and a driving term under time- 
dependent similarity transform: 



H{9) = H Q + U m (9)h m U{9) 



(1) 



where is a time-dependent rotation parameter and Hq is 
a Hamiltonian with two degenerate ground states of zero 
energy. Suppose one ground state is \ip a ) = Y^i=i a i 
and another is \ip b ) = \j m ) where {|7i) ••• l7n) > l7m)} 
spans the effective Hilbert space of interest. The 9- 
dependent term acts on a two-dimensional space spanned 
by {|7m) j It;)} f° r some 1 < j < n. That is, we define 
h m = K \^ m ) (7m | for some positive constant K and U (9) 
is a rotation in {\^ m ) , 
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Here J denotes the generator of rotation. One can check 
that 

1 ™ 
|0) = ( aj tan fl | 7m ) +^a t | 7t )), (3) 

Jl + a 2 tan 2 9 i= i 

is a ground state of H(9) with zero energy for any 9. We 
will consider 9 between and ■§. 

In adiabatic evolution, the transition rate from a 
ground state |0) to an excited state |fc) is well-known 



to be bounded by |(fc| ^§ |0) /{E k - E ) 2 \ HI- To keep 
this quantity constant as we vary 9, we write ^ = ^-^f 
and it follows that we need 
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Our model offers at least two advantages if we wish to 
adiabatically manipulate a state of the form ([3]). 



Claim 1 Let g k {9) 
9' k {0) * 0. 



\E k {9) - E {9)\. For <9 < 



Proof: If \k(9)) is an eigenstate of H{9) with energy 
E k (9), by first-order perturbation theory, 

where dH Jf^ can be explicitly evaluated as 
dH{9) 
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One readily checks that for < 9 < -| the expectation 
value of this operator vanishes if and only if the |7 m ) and 
components of the state satisfies the tan 9 ratio. If 
an eigenstate of H (9) satisfies this, it would have to be 
the ground state which stays at zero. This means the gap 
only shifts in one direction. ■ 



Claim 2 |(fc|M| )| <O(y/grf0)). 

Proof: The idea of the proof is that ^ |0) couples to a 
state with high energy, and this limits its coupling into 
low lying states during the evolution. Using ((5]), we ob- 
tain 



dH , , 
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The energy of an excited state is g k {9) = \{p\k)\ K + 
(fc| Ho \k), from which we see that |(p|fc)| < \J g k (&) j K . 
■ 

The first fact enables us to estimate efficiency by looking 
only at the end points. The second allows a higher speed 
for d9/dt when the gap is small. A more detailed estimate 
yields: 
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where g\ is the first excited state energy of Hq. We see 
that the quantities A and c m k can further suppress the 
transition rate as they may contain factors inversely pro- 
portional to the size of the system. 

To see why this is useful, consider the case 
for i = l..n and 8 running from to ? adiabatically. It 
is clear that we can turn the state -4=(|'yi) + ••• + |7n)) 
to ■^ 7 = t (|t'i) + ... + |7n> + |7m», thus adding an ele- 
ment to the uniform superposition. Reverse the process 
and we will delete the element from the superposition. 
Clearly, repeated applications of such processes with dif- 
ferent driving terms allow us to prepare a large class of 
states. Generally speaking, if we wish to create a partic- 
ular superposition of certain basis states, we can look at 
the Hamiltonian interactions we have control over that 
act either diagonally or off-diagonally on the basis. This 
set of interactions can then be thought of as a generating 
set for group elements that we want to eventually appear 
in the superposition. In the next section we will consider 
some examples to illustrate this idea. 



II. THE TORIC CODE AND THE CLUSTER 
STATE 

The ground state of the toric code Hamiltonian 0] can 
be generated efficiently by adiabatic linear interpolation 
Q. The minimal gap in the evolution is estimated to be 
of order 0(n~^). Here we show a different path where 
the gap remains constant but the evolution time is also 
0(\Zn), thus also saturating the Licb- Robinson bound 
[Ig |. Following the basis choice of the Hamiltonian 
is: 



H 



star 



z®z®z®z 



- ^2 X(g>X(g>X(g)X (7) 

plaqucttc 



where X and Z are Pauli matrices and the spin lattice 
is as shown in FigQJa). The spins live on the links (as 
this is a Z2 lattice gauge theory); a "plaquette" refers 
to the four spins around a square while a "star" means 
the four spins around a vertex. We will work within 
the ground states of the star terms so we can effectively 
ignore them. Consider the spins {1,2,3,4} in FigQJa) 
and a ^-dependent Hamiltonian: 

H p (0) = l + {sm 2 9-cos 2 6)(Z 1 +Z 2 + Z 3 + Z 4 )/4 
-2 sin9 cos 6X^2X3X4, (8) 

We can check that within the subspace of |0000) and 
1 1111) this has exactly the form of the driving term in 
© and that it takes |0000) to ^(|0000) + jllll)) as 6 
goes from to |. Other states in the Hubert space, re- 
gardless of energy, do not affect the result because they 
are decoupled from the subspace, so the effective gap is 
2 throughout the process. Since the ground state of the 




FIG. 1: (a) Interactions that appear in the toric code Hamil- 
tonian. (b) Interactions in the cluster state Hamiltonian, fol- 
lowing the same notation. 



Hamiltonian j7]) is known to be the superposition of all 
"closed strings" (plaquette excitations), the goal is to 
do the same operation of rotating the Z operators into 
XXXX for every plaquette. Consider next the plaquette 
{4, 5, 6, 7} after we have evolved {1, 2, 3, 4}. The Hamil- 
tonian we use would have to be a little different since "4" 
is already entangled: 



H' p (0) = 1 + (sin 2 9 - cos 2 0){Z 5 +Z 6 + Z 7 )/ '3 
-2 sine cos 6X 4 X 5 X 6 X 7 



(9) 



It is easy to verify that this Hamiltonian takes 
^(|0000) + |1111)) 1 000) to i(|0000000) + 10001111) + 
11111000) + 11110111)) as 6 goes from to f . Generally, 
we can add a plaquette excitation to the superposition if 
at least one out of the four spins is unentangled. Bearing 
this in mind, we can perform these adiabatic rotations on 
plaquettes in parallel (coefficients in the Hamiltonian can 
be also adjusted to give a more translationally symmet- 
rical appearance). For example, on a square lattice, we 
can operate on alternate rows of plaquettes in parallel. 
Dividing each row into odd and even plaquettes, this can 
be accomplished in two time steps. Then we operate on 
the remaining rows sequentially - to make sure we have 
an unentangled spin for every operation - which takes L 
steps if the lattice has n = L x L spins. The genus of the 
surface does not really affect this algorithm, except we 
have to notice that because the same loop can be viewed 
from two sides in a torus, there will be some plaquette 
operators that do not lead to new states in the super- 
position. As for other topologically inequivalent ground 
states, they can be generated in this manner by reversing 
the Z operators around an incontractible loop at the be- 
ginning. Compared to the linear interpolation, our evolu- 
tion path is certainly more complicated. But it provides 
an interesting perspective, as it does not rely on duality 
mapping to an Ising model. In principle the technique 
we use can be generalized to prepare any string-net con- 
densed state [1] because they are uniform superpositions 
of group elements generated by local string-nets. 



The same trick also works for the cluster state [13|. The 
cluster state used in measurement-based quantum com- 
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puting, if put on a square lattice, is the superposition 
J2 Z ={o i}(~l)'' \ z i--- z n) where r is the number of ad- 
jacent "l"'s in the configuration {zi...z n }. It is stabilized 
by the operators X) s ^s(®s' ^a') where s' denotes sites 
adjacent to site s and the tensor product is over all four 
adjacent sites. Since each term takes half the states in 
the superposition to the other half, they play a role sim- 
ilar to the plaquette operators in the toric code. We can 
set the Hamiltonian for the sites {1, 2, 3, 5, 6} in FigQJb) 
to be something similar to (JSj) : 

H c (9) = 1 + (sin 2 6 - cos 2 8)Z 3 

~2smdcosdX 3 Z 1 Z 2 Z 5 Z 6 (10) 

and take 9 from to f . The same operation can be 
repeated on sites {3, 4, 6, 7, 8}, for instance. But since we 
always have an unentangled spin every time we turn on 
a stabilizer term, unlike the toric code case there is no 
need to operate sequentially. Therefore we can operate 
on all the odd sites in the first time step and all the even 
site in the second time step, completing the algorithm in 
0(1) time. 

Note that nothing prevents us from applying pertur- 
bation gadgets 14] on the stabilizers and turning our 
Hamiltonian into a nearest-neighbor one. Thus this is 
comparable with the linear interpolation result of 
with the difference that we are preparing the cluster state 
rather than a state that can be used as one. 



IV. THE HISTORY STATE 
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FIG. 2: Comparison of energy gap between stepwise adiabatic 
rotation to linear interpolation in the generation of history 
state with L=6. For the stepwise adiabatic rotation, the gap 
never crosses the line corresponding to the final gap. (Color 
online) 



ing term in Hf\ 

H t {0 t ) = sin 2 6l t |7 t )(7 t |+cos 2 6> t |7 t+1 )(7 t+1 | 

- sin 6 t cosdt(\j t ) <7t+i| + l7*+i) (7t|)( 12 ) 



In all our examples so far the effective energy gap re- 
mains constant, so essentially we are implementing quan- 
tum circuits in the sense of [6[ . Now let us consider an ex- 
ample where the gap shrinks as the system grows. In 
it is shown that the history of a quantum circuit - a uni- 
form superposition of intermediate states, appropriately 
made orthogonal - can be generated via an adiabatic lin- 
ear interpolation. Up to small corrections unimportant 
for our purpose, the initial and final Hamiltonians are: 

L 

Hi = $> t ><7 t | 
t=i 

1 

H f = o zZ lT*> <T*I + |7t+i)<7*+i| 
z t=o 

- l7t) <7t+il - l7t+i) (7t| (hi) 

where |7t)'s represent the orthogonal intermediate states 
at time steps labeled by t. This adiabatic evolu- 
tion takes the starting state I70) to the history state 
^ *+i S{=o It*}- An obvious alternative to linear interpo- 
lation from Hi to Hf would be to apply adiabatic rotation 
step by step and rotate each term in Hi to a correspond- 



where as usual, Ot goes from to \ at each time step. 
Figf2]shows the energy gap as a function of time for L = 6 
and compares this stepwise method to the linear inter- 
polation. We can see that the minimal gap occurs at the 
very end. In 5] sophisticated techniques were used to es- 
timate the scaling of the minimal gap as a function of L. 
It is much easier in our case, even compared to the sim- 
pler analysis of [ll[. Observe that the final Hamiltonian 
is the discrete version of a kinetic energy term for a free 
particle in a box of size L. For large L, its spectrum sim- 
ply consists of standing waves with energy proportional 
to the square of momentum, which goes as the inverse 
of wavelength. Thus we immediately conclude that the 
gap is 0(L~ 2 ). Now consider the speed of evolution as 
measured by dO/dt for fixed error rate. Naively, it would 
scale as the square of the gap, i.e. 0(L~ 4 ). But we can 
use the bound in Eq[6] to get a better result. At the 
minimal gap, g\ ~ 0(L~ 2 ) and A ~ (9(L _1 ). Since the 
excited states are standing waves, it is easy to see that 
c 2 nk ~ 0(L _1 ). Altogether, we obtain dO/dt ~ 0(L- 2 ). 

Actually there is a more elegant way of using adiabatic 
rotation in this example. So far we have been taking 6 
from to ?. If we instead consider 9 evolving from J 
to j, we would take a single basis state to a uniform 
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superposition. This suggests that we try: 

H{9) = cos 2 0| 7o )(7o|+sin 2 0| 7l )( 7l | 
-sin^ cos 0(|7q) (71 1 + (71) (70 1) 
1 

+ 2 I 7 *) ^ 7 *l + l 7 * +1 ^ ( 7t+1 l 
- |7t) (7t+i| - |7t+i) (7*1 



(13) 



Here I70) plays the role of |7 m ) in Eq[3J by rotating 9 
from § to ^, we achieve the desired evolution, remark- 
ably, leaving most terms constant. The only thing we 
need to check is the energy gap at the beginning, where 
the Hamiltonian is a kinetic energy term plus a delta 
function potential of two energy units at one end. This 
is still easy to analyze compared to the linear interpola- 
tion. For sufficiently large L, the constant potential just 
amounts to a zero boundary condition for the low lying 
wavef unctions. The lowest energy state in the "box" is 
a standing wave with wavelength 4L, so the gap is also 
0(L~ 2 ). If we wish to implement the mirror evolution 
described in Q to get from the history state to |tl), 
a similar rotation can be used on the other end of the 
"box". 



V. QUANTUM SEARCH AND OPTIMIZATION 

The preparation method above generates the ground 
state of a Hamiltonian; as such it is suggestive of a 
problem-solving algorithm. Indeed, let us consider an ex- 
ample that can be analyzed exactly. In Q an adiabatic 
search algorithm was presented using linear interpola- 
tion to an oracle-like problem Hamiltonian. It turns out 
that one can formulate an alternative version using adia- 
batic rotation. The only extra resource we need is a non- 
computational state that the problem Hamiltonian anni- 
hilates (we may, for instance, add a non-computational 
state to each qubit to form a qutrit). Call this state \i) 
and our setup looks as follows: 

H{9) = (/ - |/7T,) (tti|) + sin 2 6^ |Vo> <^o| + cos 2 |i) (i| 
-sin0cos0(|V>o)(i| + N)(^o|) (14) 

where we have followed the notations of 0: \m) is 
an unknown solution state we want to obtain and \tpo) 
is a known state that has non-zero overlap with |m). 
The identity operator "I" only acts on the computa- 
tional space and is zero on \i). Denote the overlap 
by a = (xjj \m) (i.e. clq = — h= in [2j). Adiabati- 



■Sn 

to evolves the state 



to 

a 2 ,, which 
This ap- 



cally changing 9 from § 
|m). The energy gap g\{9) = 1 — sinflyT"^ 
at 6 = I is 0(a 2 ) - 0{N~ X ) if a = 
pears different from the case in [2| where the minimal 
gap is 0(N~2). We must however take into account 
|(fc| ^M- |0)|. The first excited state is proportional to 



a |m) + (sin 6*^1 - a 2 , - 1) \ip ) 
plicit calculation yields: 



. . dff . , 



< ao 



sin 9 



\/l -a, 



1 - sin 9^/1 



O(a ) 



(15) 



Now use Eqf5]to set the rotation speed: d9 ~ (gf/ao)dt. 
Integrating this from 9 = ^ to 9 — and setting 
a = — =, we obtain the total evolution time T ~ 



(VN - - l)- 1 - 0{VN) for large N. Quadratic 
speedup is achieved as promised. 

In principle we may replace / — \m) (m\ in Eq ll4l by 
other problem Hamiltonians to find other ground states. 
This would work as long as the gap at \ is sufficiently 
large. But if we want to solve NP-complete problems 
as in [l[, notice that | -0o ) (V'o I is generally nonlocal and 
does not respect the bit structure [la ]. Suppose we set 
\i) = 1 222 . . .) by using qutrits with an extra state |2) and 
assume |m) = |000...) without loss of generality, it can 
be proved that no local Hamiltonian can have a unique 
ground state of the form 



1 



= (|000. 



o tan6l |222...)) 



by arguments in Q. Since locality is important to 
any realistic algorithm, we need a local version of adi- 
abatic rotation different from (jXJ)- - For example, if 
I V'o) = ®„ 1+) where n is the number of qubits and 
|+) = ^t|(|0) + |1)), we may want to replace the driv- 
ing term by J2n ^1+) (+1 ^ where U rotates between 
|+) and 1 2). We may then expect the ground state to 
contain a piece proportional to ®„(|+) + tan 9 |2)). The 
appearance of an extra (3™ — 2") states, however, can 
significantly change the spectrum. More work would be 
required to find a feasible local model. 



VI. DISCUSSION 

The issue of locality arises in the previous section be- 
cause our rotation method requires the initial state to 
be separated from the computational states. The same 
requirement applies to the history state preparation in 
Section IV, but there locality was preserved by use of 
clock qubits (which ensures the orthgonality of the |7t)s' 
- we refer the readers to 0] for details). This suggests 
that clever use of ancilla would be important to the de- 
sign of a local adiabatic rotation algorithm. We should 
note that the form of the driving term ([2]) is largely cho- 
sen for simplicity and can be generalized. The underlying 
idea is that when we have an efficient way of adiabati- 
cally evolving from state A to state B, we may combine it 
with a time-independent "diffusion" Hamiltonian H$ to 
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attain a more complicated state B' - a state that overlaps 
B but we otherwise may not know how to get to. Such a 
perspective can be useful for algorithm design. 

Results in this article also raise another issue. In al- 
most all our examples, adiabatic rotations turn out to 
have the same efficiency as the linear interpolations - not 
only up to the same complexity class, but down to the 
order of polynomial. One might wonder if there is some 
deep connection between different adiabatic paths under 
the same set of Hamiltonian interactions. If such connec- 
tion exists, might it be possible to prove the efficiency of 
a linear interpolation path by proving the efficiency of a 
long-winded path that is easier to analyze? At this point 
this is an entirely open question. 
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